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Abstract 



We propose an approximate solution of the time-dependent Schrodinger 
equation using the method of stationary states combined with a vari- 
ational matrix method for finding the energies and eigenstates. We 
illustrate the effectiveness of the method by applying it to the time 
development of the wave-function in the quantum-mechanical version 
of the inflationary slow-roll transition. 
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There have been various approaches to the calculation of the time devel- 
opment of the Universe in the early stages of inflation, in which the inflaton 
field tp evolves from the initial unstable vacuum state in which (ip) = to 
the final stable vacuum in which (if) = ±a, say. 

At the level of quantum mechanics the earliest approach pQ used the 
Hartree-Fock method, which, although not very accurate, does give a qualita- 
tive picture of the true time development. Several years later Cheetham and 
Copeland [2] went beyond the Gaussian approximation by using an ansatz 
which included a second-order Hermite polynomial. This represented an im- 
provement on the Hartree-Fock approximation, but still did not reproduce 
the first maximum in (x 2 ) of the exact wave-function. A further attempt to 
solve this difficult non-perturbative problem involved the use of the linear 
delta expansion 3J. This was successful to the extent that it stayed close to 
the exact solution for longer than previous methods, but there seemed to be 
a barrier to its implementation beyond the first maximum. The results of 
these various approaches to the calculation of (x 2 ) are shown in Fig. 1, where 
they are compared with the exact result of Lombardo et al. obtained by 
Fourier transforming to frequency space and then back again. 
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Figure 1: (x 2 ) 1 / 2 versus t for the slow roll potential of Eq. ©. The labels HF, CC 
and LDE correspond to the Hartree-Fock method pQ, the improved Hartree-Fock 
method |2j and the fourth-order linear delta expansion [3j respectively. These are 
compared with the exact result of [I]. 

The present paper represents another approach to the problem, using 
the method of stationary states, combined with a matrix method for finding 
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the energy eigenvalues of the full Hamiltonian, together with their corre- 
sponding eigenfunctions. The method consists of converting the eigenvalue 
equation into a matrix equation by taking matrix elements with respect to 
harmonic oscillator wave functions of arbitrary frequency Q and truncat- 
ing the matrix at some finite size N. As iV increases one may expect the 
method to become more and more accurate for the lower-lying energy lev- 
els. However, for arbitrary Q the convergence is rather slow. A significant 
component of our method is to choose an optimal value of Q by a partic- 
ular version of the principle of minimal sensitivity (PMS)[5 , which greatly 
accelerates the convergence. The method has been used before for the deter- 
mination of the eigenvalues [H], but here we also make use of the corresponding 
wave-functions to implement the method of stationary states to calculate the 
time-dependence of the state from its initial configuration. It turns out that 
this simple method is extremely accurate at even quite small orders, and has 
no difficulty with the long-time behaviour. We have concentrated on this 
particular problem to exemplify the power of the method, but its range of 
application is obviously much wider. 

As stated above, we tackle the problem of solving the energy eigenvalue 
equation 

H^ n = E n ^ n (1) 

by converting it to a matrix equation, using an orthonormal basis of wave 
functions of the quantum harmonic oscillator, depending upon an arbitrary 
frequency Q = a 2 : 

<p n (x) = N n e~ a2x2 ' 2 H n (ax) , (2) 

with N n = (a/(2 n n! y/n))^. The slow-roll Hamiltonian involves a double- 
well potential, and has the specific form 

H = t^P 2 + A(x 2 -a 2 ) 2 /24 + const. 

2 1 2 2 , 4 fn\ 

= -p - -m x + gx , (3) 

with m 2 = Aa 2 /6 and g = A/24. 

Of course the infinite-dimensional matrix H n £ must be truncated to some 
finite dimension, say N x N, and then its eigenvalues can be calculated by 
simple matrix diagonalization. A brute-force approach to the problem is to 
stick with a given Q and rely on larger and larger values of iV to approach the 
desired accuracy for the lower-lying energy levels. However, much improved 
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accuracy for even modest sized matrices can be obtained by a judicious choice 
of fl The criterion we shall adopt here, which is essentially that adopted 
in [Hj, is the principle of minimal sensitivity jH] applied to the trace of the 
truncated matrix. 

The rationale behind this principle is that the eigenvalues and other exact 
quantities of the problem are independent of Q but any approximate result 
coming from the diagonalization method for finite N exhibits a spurious 
dependence on the oscillator frequency. This also applies to the trace of 
the matrix, i.e. the sum of the eigenvalues. A well-motivated criterion for 
choosing f2 is therefore to take it at a stationary point of Xjv = Y,n=o H nn , 
so that this independence is respected locally. Thus we impose the PMS 
condition 

The reason for applying this condition to the trace is that is a simple 
quantity to evaluate, and moreover it is invariant under the unitary trans- 
formation associated with a change of basis. Once Q is so determined, one 
obtains an approximation to the first N eigenvalues and eigenvectors of H by 
a numerical diagonalization of the truncated N x N matrix. One could also 
contemplate applying the PMS to the determinant, which of course shares 
the property of invariance under unitary transformations, but this would be a 
much more cumbersome calculation, and could well introduce many spurious 
PMS solutions. 

In order to implement the method we need the harmonic oscillator matrix 
elements of x p . Closed formulas have been given in Ref. jjj, which we adapt 
here for completeness. 

(x i r) _ V^! mi tC A) W (5) 

1 ' {2a) 2r ^ 2 2r ' k ~\r- X - k)\(n - k)\{2X + k)\k\ { } 

for t — n = 2A, and 



2r+1 _ V^w. 



(2a) 2r + 1 



min(n,r— A) 

E 

fe=0 



(2r+l)! 



2 2r-fc-A+±( r _ A _ k y^ n _ fc ),(2A + 1 + k)\k\ 

for t— n = 2A + 1. These formulas assume i > n, but the matrix is symmetric. 
In both cases r must be greater than A. For all other values of p, n and t the 
matrix elements vanish. 
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In addition to these we will need the matrix elements of p 2 , which are 
given by 



(P 2 



1 



-a 



^£(£ - l)8 l>n+2 - {21 + l)8 ljtl + y]{l + l){t + 2)8^ 2 



(7) 



Using these results we can write down the matrix elements of Eq. © in the 
basis of the wave functions of Eq. (j2J). The trace to order N is given by 



*T w = tf(n + =L) + £(l + 2W >), 
which turns out to have a single minimum, located at 

U PMS — -X - j^, 



where 



X = 



N 



9g(l + 2N 2 ) + V3 (N 2 m & + 27g 2 (l + 2N 



in 

2\2\ 2 



(8) 



(9) 



(10) 
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Figure 2: Trace of the truncated matrix for the Hamiltonian © normalized by the 
number of states as a function of the variational parameter 17 for a = 5, A = 0.01. 
The three curves are for N = 6, 8, 10. 



A graph of Tn/N against Q is given in Fig.Elfor various values of N. Here 
we have taken the parameters as a = 5, A = 0.01, in accordance with previous 
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work on the subject^ 12 IS]- By taking Q at the minimum, in accordance 
with PMS, we obtain a good approximation to the spectrum for even quite 
small matrices. 

For the present purposes, the important point of this method is that 
in addition to the energy levels it gives an approximation to the energy 
eigenf unctions, in the form 

N-l 

1pn(x) = d ^k ¥k{x) ■ (11) 
fc=0 

where d n k denotes the kth component of the nth eigenvector of the truncated 
Hamiltonian matrix. 

These can be used to implement the method of stationary states. Namely, 
if the initial wave-function at t = is expanded as 

oo 

*{X,0) = Y.^nix), (12) 
n=0 

its value at a later time is given by 

oo 

^(x,t) = Y a ne- lE " t Mx) (13) 

n=0 

As an intermediate step we expand ^(x, 0) in terms of the eigenvectors (p n {x)~- 

N-l 

V n {x, 0) = Cn <Pn{x) ■ (14) 
fc=0 

By comparison with Eq. (fT2"J) we see that 

N-l 

a n = ( d_1 )fc • ( 15 ) 

1=0 

Notice that by this stage Q has been determined by the PMS condition, so 
that the matrix inversion here is numerical, rather than symbolic. 

It remains to determine the coefficients c n in Eq. (|14j). The initial wave 
function used in previous studies of slow-roll inflation is given by 

H)(x,0) = ^y /4 e — 2 / 2 . 

By orthonormality, 

c n = J <f>* n (x)*{x,0)dx = N n (-) 1^ e- p2x2 / 2 H n (ax)dx, (16) 



where (3 2 = m + a 2 . By means of the change of variable y = ax 



c n = ^(-) / e-^ 2 H n (y)dy. (17) 

a V 7T / J-oo 



and the expansion of H n (y), namely 

[n/2] 



ff»G/) = EC- 1 )* ? n 9 ai^ "-V- a 

(n-2fc)!fc! 



we obtain 



where 

n-2fc+l 







J k = I y n - 2k e-Z? 1 ' dy = / z n - 2k e- z2 dz . 



oo 



-oo 



In fact c n = unless n is even n = 2£, and then 



(20) 



The coefficients a n are now given by Eq. (|T5]l . 

For comparison with previous work, we use our time- dependent wave 
function to calculate (x 2 ), given by 

(x 2 ) = J V*(x,t) x 2 V(x,t) dx (21) 

assuming that ^/(x,t) remains normalized, which we have checked. In terms 
of the ip n (x) this is 

(* 2 ) = E < a * e ~ iUJntt [ x 2 if> t (x) dx, (22) 

n,i 

where u n £ = E n — Eg. Using Eq. (|TT|) this becomes 

( x2 ) = E a n a e J2 d nkdej(x 2 )kj cos^wt . (23) 

n,£ k,j 
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Figure 3: (x 2 ) 1 / 2 versus t for the slow roll potential of Eq. Q calculated 
with our method for N=Q, 8 and 10. On the scale of the figure, the result 
for iV=10 is barely distinguishable from the exact result. 

The result for (rr 2 ) 1 / 2 is plotted in Fig. El As can be seen, the method is 
vastly superior to the previous approximative methods, and by N = 10 the 
curve can hardly be distinguished from that obtained using Fourier trans- 
form methods. Although for the present purposes a basis of size 10 was 
sufficient, considerably greater accuracy for both energies and eigenf unctions 
can readily be obtained by going to larger N. This costs little additional 
effort because, once the PMS has been applied, the Hamiltonian matrix is 
fully numerical and the calculation of its eigenvalues/eigenvectors can be cal- 
culated numerically with accuracy and speed. For example, with iV = 100 
we obtain the energy of the ground state of the quartic oscillator Hamilto- 
nian p 2 + x 2 + 2gx 4 with g = 1000 to an accuracy of 58 significant figures. 
Although we have no analytic proof, our numerical work strongly suggests 
that the error decreases exponentially with N. 

In both these cases we have been dealing with symmetric potentials, but 
with a slight modification it can easily be extended to asymmetric potentials. 
In such instead of using a basis of harmonic oscillator wave-functions 

centered on the origin, we can take them to be centered on a shifted position 
x = a. The method then contains two variational parameters Q and a, which 
are to be determined by the PMS condition that T/v be stationary in both. 

Another possible extension of the method might be to use the eigenfunc- 
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tions of the g-deformed harmonic oscillator as a basis instead of those of the 
simple harmonic oscillator. A variational method using these eigenfunctions 
has been shown to give very accurate results for the ground-state energy of 
the anharmonic oscillator [8j. 

In conclusion, we have shown that the matrix method, combined with 
the principle of minimal sensitivity applied to the trace of the Hamiltonian 
matrix, is a powerful tool for finding the spectrum of arbitrary polynomial 
potentials having only bound states. Since the method also gives the energy 
eigenfunctions it is ideal for implementing the method of stationary states 
in order to track the time development of a given initial wave function. For 
quite modest-sized matrices the method gives good accuracy, even for long 
time scales, as we have demonstrated for the slow-roll potential. 
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